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We present a study of correlations among conserved charges like baryon number, electric charge 
and strangeness in the framework of 2+1 flavor Polyakov loop extended Nambu-Jona-Lasinio model 
at vanishing chemical potentials, up to fourth order. Correlations up to second order have been 
measured in Lattice QCD which compares well with our estimates given the inherent difference in 
the pion masses in the two systems. Possible physical implications of these correlations and their 
importance in understanding the matter obtained in heavy-ion collisions are discussed. We also 
present comparison of the results with the commonly used unbound effective potential in the quark 
sector of this model. 

PACS numbers: 12.38.Aw, 12.38.Mh, 12.39.-x 



I. INTRODUCTION 

It is now well established that strongly interacting matter may exhibit a variety of phases 
depending on the ambient thermodynamic conditions. We are on the way to draw the phase 
diagram of quantum chromodynamics (QCD) which is the theory of strong interactions. The 
difficulty that we still encounter is to work with relatively strong coupling strengths in the the- 
ory. The best way to go about is to perform numerical simulation on the discretized version of 
QCD - the so called Lattice QCD (LQCD). This formulation however has not got rid of all its 
inherent technical problems and it would take some time to find the final answers fH-fl4T| . Mean- 
while one can look into the properties of strongly interacting matter through effective models 
of QCD. Polyakov loop extended Nambu-Jona-Lasinio (PNJL) model is one such model that 
successfully captures various properties of strongly interacting matter (l5l - [3Cj | . The question as 
to exactly what extent can this model emulate QCD is still a matter under investigation. The 
best way to judge it is to measure several sensitive quantities in this model and contrast some of 
them to that available in the LQCD measurements. The correlations among conserved charges 
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are some such quantities that we intend to investigate here. Essentially, the fluctuations and 
correlations of conserved charges and their higher order cumulants provide information about 
the degrees of freedom of strongly interacting matter. These can be extracted from the PNJL 
model and LQCD through the study of diagonal and off-diagonal susceptibilities respectively. 
They can also provide information about the existence of critical behavior, if any. This phe- 
nomenological study is important as there is a rapid progress in the experimental front at the 
facilities at CERN and Brookhaven heavy-ion colliders. 

On the Lattice, some of these correlators have been measured at zero chemical potential. 
With 2 flavors it was shown [3lT[34| that the fluctuations rise rapidly around the crossover 
region from hadronic matter to quasi-free quark matter. The higher order cumulants show 
non-monotonic behavior 0, Hi| . Similar measurements have been made in 2 flavor PNJL model 
with three-momentum cutoff regularization [U [3tl - [38j . Quark number susceptibility (QNS) 
at finite density has been estimated in some works within 2 flavor PNJL model [39|]. Also, 
QNS has been studied in Hard Thermal Loop approximation [40M42| . Recently the idea of the 
numerical Taylor expansion in terms of chemical potential for PNJL model has been used within 
the constraint that the net strange quark density is zero, which is the case in ultra relativistic 
heavy ion collision [43j . 

For 2+1 flavors, fluctuations have been recently measured in LQCD 
PNJL model both with the usual unbound effective potential (UEP) [l8|, 
bound effective potential (BEP) [49|. Similar calculations have been carried out in Polyakov 
loop coupled quark-meson (PQM) model (50l - [53j and its renormalization group improved version 

m 

In this work we investigate the off-diagonal susceptibilities which give the correlations among 
different conserved charges. Our paper is organized as follows. In Sec. HH we discuss the 
basic formalism of the PNJL model as well as the method of extracting the Taylor expansion 
coefficients of pressure that gives the various susceptibilities. In section IIIII we present and 
discuss our results together with a comparison with the data obtained in LQCD. The last 
section contains a summary and our conclusions. 
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II. FORMALISM 



A. Thermodynamic Potential 

2+1 flavor PNJL model with unbound effective potential has been studied elaborately in 
a number of recent works [HI, [l|| [2(| [23j]. To introduce a bound in the effective potential 
eight quark interaction terms have been introduced in 2+1 flavor NJL model |2J-|27[ and in 
the 2 flavor PNJL model [H, |2{|. We developed the 2+1 flavor PNJL model with bound 
effective potential to study finite temperature and chemical potential properties (30j within 
three-momentum cutoff regularization scheme. Comparing with available LQCD data the model 
was shown to reproduce various aspects of QCD thermodynamics quite satisfactorily. We shall 
be using this model in the present work. The relevant thermodynamic potential in the mean 
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field approximation can be written as [30], 
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where <?s and c/d are the four quark and six quark coupling constant and <?i and 32 are the 
eight quark coupling constant. Here af = (ipfipf) denotes chiral condensate of the quark with 

flavor / and £7/ = wp 2 + M? is the single quasi-particle energy. Here, constituent mass Mf of 
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flavor / is given by the self-consistent gap equation; 

M f = m f ~ 2 9s°f + ^-a f+ ia f+2 - 2g 1 o i {a 2 u +a 2 d + a 2 s ) - Ag 2 a) 

where /, / + 1 and / + 2 take the labels of flavor u, d and s in cyclic order. In the above 
expression, the vacuum part integral has a ultraviolet cutoff A. For fixing the parameters m s , 
A, gs, go-, gi, gi we have used the following physical conditions (30j ; 

= 138 MeV m K = 494 MeV m v = 480 MeV m v , = 957 McV 
U = 93 MeV f K = 117 MeV 

and m u is kept fixed at 5.5 MeV. The parameters are given in table U for UEP and BEP. 
The Polyakov loop $ and its charge conjugate 5> are defined as, 

$ = (Tr c L)/N c , $ = (Tr c Lt)/N c 

where, L is the Wilson line given by, 



"Pexp i / A^dr 
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The Polyakov loop potential W with the Vandermonde (VdM) term can be expressed as [2! 

U'(<f>, §, T)/T 4 = U(<f>, §, T)/T 4 - k ln[ J($, $)] 
where T) is the Landau- Ginsburg type potential given by fl(| . 

&4 
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with, 



6a(T) = a + J) + a 2 (^) 2 + a 3 (^) 3 , 
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Interaction 
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gi x 10" 21 
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(MeV) 


UEP 


5.5 


134.758 


631.357 


3.664 


74.636 


0.0 


0.0 


0.13 


181 


BEP 


5.5 


183.468 


637.720 


2.914 


75.968 


2.193 


-5.890 


0.06 


169 



TABLE I: Parameters and Tc for UEP and BEP type Lagrangians. 



and 63, 64 are constants. To is the deconfinement temperature in a pure gauge theory. The 
VdM determinant in eqn. @ is given by (22j, 

J[$, $] = (27/24vr 2 )(l - 6$$ + 4($ 3 + $ 3 ) - 3($$) 2 ) 

Here n is a phenomenological constant which is determined by reproducing the pressure calcu- 
lated on Lattice. For the Polyakov loop potential we choose the parameters which reproduce the 
Lattice data of pure gauge thermodynamics According to pure SU(3) lattice gauge theory 
value of To is found to be 270 MeV. However we took To as 190 MeV to get the crossover tem- 
perature (T c ) consistent with the LQCD data. Various thermodynamic quantities like scaled 
pressure, entropy and energy density are reproduced extremely well in Polyakov loop model 
using the ansatz ([3]) and Q with parameters summarized below, 

a = 6.75, ai = -1.95, a 2 = 2.625, a 3 = -7.44, 63 = 0.75, b 4 = 7.5, T = 190MeV 

There have been some work [55| where the Polyakov loop parameters have been obtained by 
fitting with the full LQCD results rather than with the pure gauge theory results as done here. 
However, since we have fitted the pressure obtained in our model with that obtained in the 
full LQCD to obtain the parameter k, the effect of full QCD is incorporated. In this work the 
Polyakov loop is a global object. One can improve on it by taking in to the consideration of 
quantum and local corrections as done in [56l.l57l|. However, in this work, we are mainly looking 
at the trends of different observables rather than exact matching with the LQCD results. As 
more and more refined LQCD results are coming up we understand that exact quantitative 
status is going to change. Hence we do not incorporate those involved calculations in this work. 



B. Taylor expansion of pressure 



The pressure of the strongly interacting matter can be written as, 

P(T, fi B , fi Q , us) = -to(T, fi B , Mo, Ms), (5) 

where T is the temperature, [is is the baryon (B) chemical potential, fiQ is the charge (Q) 
chemical potential and fi$ is the strangeness (S) chemical potential. From the usual thermo- 
dynamic relations the first derivative of pressure with respect to quark chemical potential fi q is 
the quark number density and the second derivative corresponds to the QNS. 

Our first job is to minimize the thermodynamic potential numerically with respect to the 
fields a u , (id, c s , $ and $. Using these values of the fields we get the mean field value for 
pressure using the equation ([5]). The scaled pressure obtained in a given range of chemical 
potential at a particular temperature can be expressed in a Taylor series as, 

P (T,^ Q ,, S ) = ^ c ff/ { T){^n^{^f (6) 

n— i-\-j-\-k 
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where, 



c ^< fe 1 > i\j\k\ d(t^y d(^y d(^) k M „ Q , S=0 1 1 

The flavor chemical potentials [i u , /j^, fi s are related to /j^, /jq, /_ts by, 

12 11 11 

M« = 3 Mb + ttMq, = 3 Ms ~ ttMq, Ms = gMs - 77MQ - Ms (8) 

Here the odd terms vanish due to CP symmetry at vanishing chemical potential and the cor- 
relation functions with i + j + k even are nonzero. In this work we evaluate the correlation 
coefficients up to fourth order which are generically given by; 

xy = 1 d^(P/T*) 

where, X and Y each stands for B, Q and S with I / F. To extract the Taylor coefficients, 
first the pressure is obtained as a function of different combinations of chemical potentials 
for each value of T and fitted to a polynomial about zero chemical potential using the gnu- 
plot fit program [58[. Stability of the fit has been checked by varying the ranges of fit and 
simultaneously keeping the values of least squares to 10 -10 or even less. 



III. RESULTS AND DISCUSSION 

We now set out to present the results obtained for correlation among different conserved 
charges. First the leading order correlations are shown and compared with those of LQCD. 
Later we discuss the behavior of some higher order correlations predicted from PNJL model. 



Let us consider the baryon-strangeness (BS) correlation. In Fig 1(a) leading order BS corre- 
lation is shown and compared with LQCD data. Though there is good qualitative agreement 
the results have significant difference quantitatively. This is not very surprising as the inherent 
physical masses of the constituents in our model and those on the Lattice are substantially 
different - the ratio of pion to kaon mass is about a factor of 2 larger on the Lattice than the 
physical value. As we shall see this kind of departure remains for almost all the correlations 
measured on the Lattice and in PNJL model. 

The BS correlation normalized to the strangeness and baryon number fluctuations respec- 
tively are given by 

XBS _ 1 eff 
^BS — — —7: — 5- 

Xss 2 cf 
XBS 1 cf s 



C 



SB 



Xbb 2 a 



-11 

B 
2 



d 2 P d' 2 P 

where wc have used the notation; xxy = 7; ^ — and xxx = 9 . It was argued in Ref. 59j 

o fix o fly OfJ-x 
that Cbs has entirely different behavior in hadron gas and in QCP, and therefore this can 
be a reasonable diagnostic tool for identifying the nature of the matter formed in heavy-ion 
collisions through event-by-cvent fluctuations. In quark phase, baryon number and strangeness 
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FIG. 1: Leading order Baryon-strangeness (BS)correlation as a function of T/Tc- Lattice data taken 
from Ref.[45||. Arrows on the right indicate the corresponding SB limit. 



are strongly correlated through the strange quark indicating Cbs should approach its Stcfan- 
Boltzmann (SB) limit as soon as quark quasi-particles are dominant. The corresponding value 

-it 



for the ratio — 5- is -1. From Fig 1(b) we see that closely above Tq the ratio reaches its SB 



limit. 

At low temperatures in the hadronic phase the situation is different. Numerator of Cbs has 
contributions from strange baryons only, whereas the denominator has contributions from all 
strange hadrons. So the ratio approaches zero as the temperature is decreased. 

A similar behavior is found in the ratio Csb which gives the BS correlation normalized to 
the fluctuation of baryon number. This is shown in Fig 1(c) Here again the high temperature 



behavior is consistent with a quark quasi-particle picture and at the low temperatures the 
strange baryon correlation is much smaller than the baryon fluctuation due to the large mass 
of the strange baryons. For both these ratios we see a nice qualitative agreement with lattice 
data though quantitative disagreement persist. 



Going a step further we show the behavior of some fourth order correlations 



„BS „BS 



-22 



-31 



and 

cf 3 s in Figl^J At low T, in the hadronic phase, all three correlations go to zero. On the other 
hand the correlations approach their SB limit at temperature close to 2.5T C . For cf 2 and 
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FIG. 2: Fourth order baryon-strange correlation coefficients as 
right indicate the corresponding SB limit. 



function of T/Tq- Arrows on the 



c^f there are two cusps, one corresponding to the chiral transition in the light quark sector 
and the other corresponding to the strange sector. Similar features for diagonal correlators 
were discussed by us in Ref. [49j| . For cf 3 s the strange sector completely overwhelms the light 
quark sector as expected, and so there is only one peak at 1.5T C . Therefore, if somehow these 
three correlations freeze out earlier than thermal and chemical freeze-out in heavy-ion collision 
experiments, they will not only be good indicators of the crossover but can also draw the explicit 
distinction between the chiral transitions in the light and strange quark sectors. 

We now turn to baryon-charge (BQ) correlation. In Fig 3(a) leading order correlation is 



shown. At both very low and at very high temperatures c n is zero. This is because at low 
temperatures the contributions from heavy baryons decrease. On the other hand in the high 
temperature weakly interacting phase the baryon and charge quantum numbers are completely 
independent of each other. The peak in cf± occurs slightly above T c and is a clean indicator 
of the crossover. We show normalized with respect to cf and c% respectively in other two 
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FIG. 3: Leading order baryon-charge (BQ) correlation as a function of T /Tc- 
Ref.fH. 



Lattice data taken from 



panels of FigJJl These can be expressed in terms of the following: 
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Both the ratios go to zero at high T. At low T the ratios show different behavior. While 
becomes small at low temperatures due to heavy baryons, is not so small due to the contribu- 
tions from light mesons. Thus Cqb remains non-zero whereas Cbq goes to zero. Interestingly, in 
case of a complete thermal equilibrium, in the fireball created in heavy-ion collisions, these two 
quantities, if measured, would give valuable insight into the quantitative aspects of the PNJL 
model at low temperatures in the hadronic phase. The difference in Lattice and PNJL model 
studies which is more pronounced below T c , would then have to confront this measurement in 
the experiment. 

In Fig0]we have plotted the 3 fourth order correlation coefficients c^ 2 , of® and cf 3 . All of 
them show a pronounced peak close to T c . Among these c 31 is non-zero only at T c and so is 




FIG. 4: Fourth order baryon-charge correlation coefficients as a function of T/Tc- Arrows on the right 
indicate the corresponding SB limit. 



the best indicator of the crossover. These correlators do not show the double peak structure as 
contribution from strange sector is subdominant. 

In FigJ5]the leading order QS correlation is shown. As in the case of baryon number, charge is 
also strongly correlated to strangeness through strange quarks and therefore at high temperature 
cff approaches its SB limit. At low temperatures, the strangeness carriers become heavier and 
hence the correlation decreases. Interesting features arise in the ratios, 

1 QS 

p XQS 1 Cii 

Xss 2 cf 

n _ XQS _ 1 Cti 

XQQ 2 C Q 

The ratio Csq decreases with decreasing temperature until it reaches T c . There it shows a 
plateau both in PNJL and Lattice measurements as most likely this is where hadronic degrees 
of freedom start to become dominant. Thereafter the ratio again decreases. This seems to 
indicate that as we raise the temperature, the strangeness content and hence its correlation 
keeps on increasing reaching a saturation. This is suddenly broken by the liberation of new 
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FIG. 5: Leading order charge-strange correlation as a function of T/Tq- Lattice data taken from 
Ref . [45|[ . Arrows on the right indicate the corresponding SB limit. 



degrees of freedom in terms of quark quasi-particles just above T c . More detailed investigation 
is necessary to confirm this picture. In retrospect we do expect and sec a similar behavior for 
Csb and Cqb though not highly prominent. 

At large temperatures Cqs takes up non-zero value [60l ] . However it increases with decreasing 
temperature due to the same reason as Cqb, i.e., the strangeness fluctuations decrease at a faster 
rate at lower temperatures. The saturation effect near T c is normalized out. 

In FigJS] we have plotted the fourth order correlations. All three fourth order correlations 
c 22 > c 3i an d C ?^ S show similar behavior as for the BS correlations. Therefore these two sets 
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can be used complementarily to understand the state of affairs in heavy ion collisions. 



IV. CONCLUSION 



We have studied the correlations between different conserved charges in the PNJL model. 
The baryon-strange (BS), baryon-chargc (BQ) and the charge-strange (QS) correlations were 
obtained by fitting the pressure in a Taylor series expansion around vanishing chemical potential. 
The ratio of these leading order correlators to the respective quadratic diagonal correlators were 
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FIG. 6: Fourth order charge-strange correlation coefficients as a function of T/Tc- Arrows on the 
right indicate the corresponding SB limit. 



obtained. The results were shown both for the physically bound effective PNJL model (BEP) 
as well as the conventional PNJL model without a physical bound (UEP). Lattice QCD data 
exists for the leading order correlators and were contrasted against our calculations. 

In general the comparison with Lattice data gave excellent qualitative agreement reproducing 
all the physical features in all the coefficients that could be compared. However quantitative 
differences remain which we believe is mainly due to the difference in effective masses. Interest- 
ingly the UEP seems to be closer to the Lattice results. One can therefore say that the PNJL 
model is standing quite strong as an effective model of QCD. 

As we discussed in the main text the combined studies of various correlators would give us 
a clearer picture of the matter created in heavy-ion collision experiments. The leading order 
coefficients can be most useful in identifying if the QGP is formed, while the higher order 
coefficients could identify the crossover region. 

We have noted a slight saturation of charge and strangeness through the saturation of the 
ratios Csq, Csb and Cqb in a small temperature region just below T c . This is found both in 
the PNJL model and in LQCD. Perhaps this is when the degrees of freedom crossover from 
hadronic to the partonic ones. 

The higher order correlators containing strangeness show two cusps at around T c and around 
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1.5T C , corresponding to the temperatures where the chiral crossover in the light and strange 
quark sectors occur. This widens the temperature range where the fluctuations remain much 
larger than the values at low temperature thus increasing the scope of identifying the possi- 
ble existence of the high temperature phase in heavy-ion collisions. On the other hand the 
higher order BQ correlators have a sharp peak and is an accurate tool to decide the crossover 
temperature. 

Even if we are not lucky enough to identify the high temperature phase from the correlations 
which need to get frozen much before thermal or chemical freeze-out, the various equilibrium 
thermodynamic measurements of the correlators would help us in determining the finite tem- 
perature behavior of the hadronic sector studied theoretically using PNJL model and Lattice 
QCD. 
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Note Added : After finishing this work we learned of a recent work by W. Fu and Y. Wu 
(arxiv : 1010.0892 [hep-ph]). However they have calculated the off-diagonal susceptibilities only 
for the unbound effective potential of the PNJL model. Furthermore, they have not included 
the Vandermonde term in the Polyakov loop potential. 
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